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ABSTRACT 

We report on the detailed study of the 2008 October outburst from the anomalous X- 
ray pulsar (AXP) IE 1547.0-5408 discovered through the Swift/Bmst Alert Telescope (BAT) 
detection of SGR-like short X-ray bursts on 2008 October 3. The Swift/X-ray Telescope 
(XRT) started observing the source after less than 100 s since the BAT trigger, when the flux 
(~6 x 10~ n erg cm~ 2 s _1 in the 2-10 keV range) was >50 times higher than its quiescent 
level. Swift monitored the outbursting activity of IE 1547.0-5408 on a daily basis for approx- 
imately three weeks. This strategy allowed us to find a phase-coherent solution for the source 
pulsations after the burst, which, besides P and P, requires a positive P term (spin-down 
increase). The time evolution of the pulse shape is complex and variable, with the pulsed frac- 
tion increasing from 20% to 50% within the Swift observational window. The XRT spectra can 
be fitted well by means of a single component, either a power-law (PL) or a blackbody (BB). 
During the very initial phases of the outburst the spectrum is hard, with a PL photon index 
T — 2 (or kT ~ 1.4 keV) which steepens to T ~4 (or kT - 0.8 keV) within one day from the 
BAT trigger, though the two components are likely present simultaneously during the first day 
spectra. An INTEGRAL observation carried out five days after the trigger provided an upper 
limit of -2 x lO" 11 erg cm" 2 s" 1 to the emission of IE 1547.0-5408 in the 18-60 keV band. 

Key words: stars: neutron - X-rays: bursts - X-rays: individual: IE 1547.0-5408. 



1 INTRODUCTION 

Magnetars are a small class of isolated X-ray pulsars, believed 
to have the highest magnetic fields known to date (~10 14 - 
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10 15 G; [Duncan & Thompson! fl992l : iThompson & Duncanlfl993h . 
This class comprises the Anomalous X-ray Pulsars (AXPs) and the 
Soft Gamma-ray Re peaters (SGRs), o bservationally very similar in 
many respects (see iMeregh etti 2008 for a recent review): a spin 
period in the 2-12 s range, large period derivatives (10 -13 -10~ 10 
6), unpredictable bursting activity on different timescales (from 
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ms to hundreds of seconds) and luminosities (10 38 -10 46 erg s _1 ). 
Until not long ago AXPs were thought as persistent and stable 
X-ray sources. Only in 2003 the first transient AXP was discov- 
ered by RXTE, namely XTE Jl 8 10—197, which displayed a factor 
of >100 flux enhancement with respect to the pre-outburst lumi- 
nosity l evel as seen by the ROSAT and Einstein missions (^10 33 
erg s-Mlbrahim et amooillsrael et alj2004lGotthelf et al.l2004l ; 
iGotthelf & HalperrJl2005l : lBernardini et al.ll2009h . 

Even more surprising was the discovery of a highly variable 
pulsed radio emission which followed the XTE J 1 8 1 0- 1 97 outburst 
i lCamilo et ail 120061 : iHelfand et aUl2007h. n ever observed before 
in any other magnetar l Burgav et alj 20061) . Moreover, the tran- 
sient nature of this AXP provided the first hint that a relatively 
large number of members of this class has not been discovered yet, 
and suggested that others would manifest themselves in the future 
through a phenomenology (outburst) similar to that displayed by 
XTEJ1810-197. 

Indeed, after this first discovery, two more transient AXPs 
have been detected showing similar outbursts, CXOU J164710.2- 



455216, dMuno et al.1 120071 ; llsrael et all l2007j) and IE 1547.0- 
5408. The latter AXP was first suggested as a candidate magnetar 
in the supernova remnant G327 .24-0.13 through X-ray observa- 
tions ( Gelfand & Gaensler 2007), and subsequently recognized as 
a radio transient magnetar th rough the discovery of radio pulsations 
dCamilo et aT1l2007bl. 120081) a t a period of ~2. 1 s, and of an X-ray 
outburst dHalpern et alj|2008h . In particular, Swift/X-ray Telescope 
(XRT) observations taken in Summer 2007 revealed IE 1547.0- 
5408 at an X-ray flux level of ~5 x 10 -12 erg cm -2 s _1 (more than 
one order of magnitude brighter than in quiescence), with charac- 
teristics similar to those of XTE J1810-197 in outburst. No X-ray 
bursts were observed during the 2007 outburst of IE 1547.0-5408, 
possibly due to a sparse X-ray coverage. More generally, the first 
phases of the 2003 and 2007 X-ray outbursts from XTE J 1 8 1 0- 1 97 
and IE 1547.0-5408, respectively, were missed. 

Relatively deep XMM-Newton pointings were obtained in Au- 
gust 2006 and 2007, the first during quiescence, and the second dur- 
ing the outburst decay. Pulsed fractions of about ~15% and ~7% 
were inferred for th e 2006 and 2007 observations, respectively 
dHalpern et al]|2008h . The 2006 quiescent spectrum (~4 x 10~ 13 
erg cm~ 2 s _1 in the 1-10 keV range) can be fitted well with the 
usual thermal (blackbody) plus non-thermal (power -law) model 
(kT = 0.40 keV and V = 3.2; lHalpern et alj|2008l) . with a ID 
resonant cyclotron scattering model (with a surface temperature 
fcT=0. 33 keV, an elect ron velocity /3=0.32 and an optical depth 
r=1.0; lReaetal.ll2008h . or with 3D Monte Carlo model (with a 
surface temperature of fcT=0.33 keV, an ele ctron bulk veloc ity of 
/3 butk =0A5, and a twist angle of Ad>=l. 14: IZane et alj|2009l) . On 
the other hand, the 2007 spectrum (~3 x 10~ 12 erg cm~ 2 s _1 in 
the 1-10 keV range) is characterized by a harder emission with 
kT = 0.52 keV and T = 1.8 (by using the BB plus PL model). 

Deep infrared observations taken from ESO-VLT during the 
2007 outburst have revealed four objects consistent with the ra- 
dio AXP position, although none of them showed variability 
(Mignani et al. 2009). Recently, during the giant outburst detected 
from IE 1547.0-5408 on 2009 January 22, a relatively bright tran- 
sient IR source (K s ~ 18.5 mag) was disc overed within the radio 
positional uncertainty region of the source dlsrael et alj|2009h . 

Within the magnetar scenario, AXP outbursts are thought to 
be caused by large scale rearrangements of the external magnetic 
field, either accompanied or triggered by fracturing of the neutron- 
star crust. These events may result in renewed magnetospheric ac- 
tivity through the interaction between thermal surface photons and 



charges flowing along the (closed) field lines. In addition, hot spots 
on the neutron-star surface may appear where the currents impact 
on the star. Repeated resonant cyclotron scatterings onto the mag- 
netospheric particles result in a modification of the seed thermal 
spectrum with the app earance of a non-thermal, high-energy tail 
dThompson et alj|2002h . 

During the first month since the onset of an outburst, the 
source experience a large flux variations (approximately by one 
or two orders of magnitudes). In this respect, monitoring cam- 
paigns, obtained just after the outburst onset, are giving an impor- 
tant opportunity to track the evolution of the main physical param- 
eters and, therefore, to check the goodness of the (empirical and/or 
more physical) components used to model the spectra. Here we re- 
port on the discovery of bursting activity and a new outburst from 
IE 1547.0-5408, as observed by Swift (©, which was promptly 
triggered by the onset of the bursting activity of this transient AXP 
on 2008 October 3. The Swift monitoring of the source in the 100 s- 
22 days time interval allowed us to carry out a detailed timing and 
spectral study, the first ever for this source, of the initial phases of 
its outburst on a daily basis. We also report on the results of an IN- 
TEGRAL Target of Opportunity observation performed during the 
X-ray outburst (|}3), and discuss our findings in comparison with 
other magnetar outbursts (ij4). 



2 SWIFT OBSERVATIONS AND DATA ANALYSIS 



The Swift satellite dGehrels et al.l I2004T) is an efficient observa- 
tory for the discovery and the multiwavelength monitoring of AXP 
bursting activity. Whenever the wide-field coded- aperture mask 
Burst Alert Telescope (BAT; iBarthelmv et al. 1 12005b triggers on a 
burst from an interesting hard X-ray transient, the system on-board 
can rapidly re-point the satellite in order to bring the source in 
the narrower field of view of the other two tel escopes on-board, 
the X-Ray Telescope (XRT; lBurrows et alj|2005l) a nd the Ultravio- 
let/Optical Telescope (UVOT; Roming et alj 20051) . In this section 
we present the results obtained from our analysis of the Swift BAT 
and XRT observations of AXP IE 1547.0-5408 performed during 
its bursting activity on 2008 October, being the source too absorbed 
to be detected with the UVOT. 



2.1 Burst Alert Telescope data 

On 2008 October 3 at 09:28:08 UT BAT triggered on and local- 
ized a short burst from a position consistent with that of the AXP 
IE 1547.0-5408 (trigger 330353; iKrimm et alj |2~008a) which was 
followed by several other bursts extending to at least 09:38:24 UT. 
About 2 hours later, at 11:16:13 UT, BAT detected one last, bright 
burst which reached a peak count r ate of ~20 000 count s s _1 in the 
15-350 keV band (trigger 330367; IKrimm et al.ll2008bh . The BAT 
time coverage of the source was discontinuous along the bursting 
period, totalling an effective exposure time of 4 ks (see FigureQJ. 

Datasets containing information for individual photons (event 
data) and data composed of 80-channel histograms with a typical 
integration time of ~5 minutes (survey data) were distributed. In 
the following, we present the results of our analysis for the bursts 
detected in the event data. 

The data reduction was performed using version 2.9 of the 
standard Swift/BAT software distributed within FTOOLS under the 
HEASOFT package (version 6.5.1). A weight map, taking into ac- 
count the contribution from the source to the counts in each BAT 
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Figure 1. The light curve of IE 1547.0-5408 generated on-board by the Swift Burst Alert Telescope (rate data). Black correspond to temporal ranges for which 
only rate data were available, red and blue to the temporal ranges for which also event and survey data, respectively, were distributed. To corresponds to the 
first BAT trigger (trigger 330353; 2008-10-03 09:28:08 UT). The units of the horizontal axis variable are seconds. 
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Figure 2. Swift/BXI background subtracted light curves of the IE 1547.0-5408 bursts. Time bins are of 5 ms. The units of the 
horizontal axis variable are seconds. Negative times marks the not triggered burst few seconds before To. 



detector (mask-weighting technique)Q was applied in order to ob- 
tain the background-subtracted counts of the source. 

In Figure E] the 15-50 keV light curves of bursts with more 
than 350 counts are plotted. For each burst we extracted a spec- 
trum. All spectra are fitted well by a single blackbody model (BB) 
with a temperature of about 1 1 keV and an emitting radius of about 
1 km (see Table Q}. The limited statistics of BAT spectra did not 
allow us to check whether an additional B B component is present 
similar to the case of SGR 1900+14 dlsrael et aljl2008bh . We also 
tried to fit the 15-100 keV spectra with a power law and an opti- 
cally thin thermal bremsstrahlung obtaining worse reduced \ 2 val - 

1 See 

http://heasarc.nasa.gov/docs/swift/analysis/threads/batfluxunitsthread.html 



ues (in the 1.4-M.6 and 1.7-^2.9 range for the faintest and brightest 
bursts, respectively, and for 38 degree of freedom). The characteris- 
tic time-scales r of each burst decay was also inferred by fitting the 
corresponding light curves with one or two (for the longest events) 
exponential laws and are reported in TableQ] 

Whenever the bursts had sufficient statistics, time-resolved 
spectra were generated corresponding to the rise, peak, and de- 
cay phases of the bursts. The results of the fitting procedure is 
reported in TableQ] No significant variations were found in the 
spectral parameters, except for the normalization factor due to the 
luminosity evolution (decay) during the burst. On the other hand, 
time-resolved spectra show a moderate spectral evolution with the 
peak BB temperature being systematically softer than that dur- 
ing decay (similar results where obtained for SGR 1806-20 and 
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Table 1. Spectral fit results (in the 15-100 keV range) of the bright bursts detected by Swift/BXT for the BB model. We assume a 
distance to the source of 4 kpc. To corresponds to the first BAT trigger (2008-10-03 09:28:08) and errors are reported at Icr confidence 
level. 



Burst time (T-To) r Net counts Duration kT R L Fluence Xr(d.o.f) 

(s) (s) (keV) (km) ergs - 1 ergcm~ 2 



-2.45 °- 007 lo:o02 358 0016 lal ti:o L5 ±0.3 

o 0.0181°;!^ 1492 0070 ia9 ^o:8 L3 io:2 

111.6 0.0361°;°°^ 2964 0.098 12.4i°j 1.2±°^ 

282.45 0.015+o°°i 1451 0.035 11.3+.0® 1 - 6 -ol 

427.9 o.o38i°:°o? 988 0100 9 - 3 ±ii l - 2 -a% 

524.2 °- 038 to:oo5 8081 °- 254 n - 8 io.3 lA -o\ 



0.133 



+0:009 

-0.008 



524.2 rise 310 0.010 lUllo 1 - 5 io'i 

524.2 peak 742 0.007 8.3±^ 4.6+^1 

524.2 decay 7029 0.237 H- 9 ^o3 l3 -o\ 

612.55 0.0121°;°°^ 757 0.100 H.ltJ;? 0.7t°j 

6485 0.0981°'°°^ 14317 0.730 10.8t°3 l3 -o\ 

3+0:019 

6485 rise 676 0.020 9.6±J ° 2.0±°* 

6485 peak 806 0.010 7.4±°-j 5.0l° j 

6485 decay 12835 0.700 10.9+° j 
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SGR 1900+14: IgoTz et al.ll2004l200f3 ; |lsrael etai]|2007bl) . 

The presence of persistent emission during the non-bursting 
time intervals was also searched for in the event data. An image 
was generated from all the events of each file covering time inter- 
vals where no bursts were present and the standard source detection 
was run (BATCELLDETECT). Two time intervals were investigated, 
the first one ranging from t = —41.5 s to t = 965.5 s (net exposure 
time 1 001 s) and the second one from t = 6 248.5 s to t = 7300 
s (net exposure time 1 050 s). No significant emission in the di- 
rection of AXP IE 1547.0-5408 were found. The 3a upper limits 
on the flux in the 14-50 keV band for the two intervals described 
above are 6 x 10" 10 and 2 x 10" 10 erg cm s , respectively. 



2.2 X-Ray Telescope data 

Swift executed a prompt slew on IE 1547.0-5408 and the first XRT 
observation started only 99 s after the BAT trigger. Table |2]reports 
the log of the observations that were used for this work. The XRT 
uses a CCD detector sensitive to photons with energies between 
0.2 and 10 keV. We considered Windowed Timing (WT) and Pho- 
ton counting (PC) mode data. In PC mode the entire CCD is read 
every 2.507 s, while in WT mode only the central 200 columns are 
read and only one-dimensional i maging is prese rved, achieving a 
time resolution of 1.766 ms (see lHilletai] |2004 for a detailed de- 
scription of XRT modes). 

The data were processed with standard procedures using the 
FTOOLS task XRTPIPELINE (version 0.12.0) and events with grades 
0-1 2 and 0-2 were selec ted for the PC and WT data, respectively 
(see Burrows et "al]|2005l) . For the timing and spectral analysis, we 
extracted the PC source events in the 0.5-10 keV range within a cir- 
cle with a radius of 30 pixels (~7l"). The WT data (in the 0.5-10 
keV range) were extracted in a rectangular region 40 pixels long 
(and 20 pixels wide) along the image strip. To estimate the back- 
ground, we extracted PC events within an annular region (radii of 
50 and 75 pixels) centered on the source and WT events within a 




Time (d - 54742.0 MJD) 

Figure 3. 0.5-10 keV SwiftKKF pulse phase evolution with time, together 
with the time residuals (central panel) with respect to the two phase- 
coherent timing solutions discussed in the text: P— P (triangles) and P—P- 
P (squares). The solid lines in the upper panel represents the latter timing 
solution. The line in the central panel is the cubic fit to the P-P solution. 
The pulsed fraction (defined as the semi-amplitude of modulation divided 
by the mean source count rate) evolution is reported in the lower panel 
(where 0.2 means 20% pulsed fraction). 



rectangular box (40 x 20 pixels) far from IE 1547.0-5408. 
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Table 2. Swift/XRS observation log of IE 1547.0-5408. 



Sequence Mode Stall time (UT) End time (UT) Exposure 51 Time since trigger 

yyyy-mm-dd hh:mm:ss yyyy-mm-dd hh:mm:ss (ks) (d) 



00330353000 b 


PC 


2008 


10-03 09:29:47 


2008- 


10-03 12:52:41 


3.9 


99 s 


00330353001 


WT 


2008 


10-03 12:54:35 


2008- 


10-04 01:39:47 


14.2 


0.14 


00330353002 


WT 


2008 


10-04 17:27:08 


2008- 


10-04 20:53:52 


4.8 


1.33 


00330353004 


WT 


2008 


10-05 09:42:43 


2008- 


10-06 21:10:00 


10.5 


2.01 


00330353005 


WT 


2008 


10-07 06:39:12 


2008- 


10-07 18:15:00 


7.6 


3.88 


00330353006 


WT 


2008 


10-08 00:20:12 


2008- 


10-08 05:13:00 


4.5 


4.62 


00330353007 


WT 


2008 


10-09 05:20:33 


2008- 


10-09 13:30:31 


1.0 


5.83 


00330353008 


WT 


2008 


10-10 07:05:50 


2008- 


10-10 15:16:00 


3.9 


6.90 


00330353010 


WT 


2008 


10-12 07:11:17 


2008- 


10-12 18:31:00 


3.7 


8.90 


0033035301 1 


WT 


2008 


10-13 04:05:15 


2008- 


10-13 15:23:00 


3.4 


9.77 


00330353012 


WT 


2008 


10-16 01:05:20 


2008- 


10-16 06:09:00 


4.0 


12.65 


00330353013 


WT 


2008 


10-17 23:38:20 


2008- 


10-18 14:13:00 


3.7 


14.59 


00330353014 


WT 


2008 


10-20 14:47:50 


2008- 


10-21 00:00:00 


3.9 


17.22 


00330353015 


WT 


2008 


10-22 08:03:14 


2008- 


10-22 13:08:00 


3.9 


18.94 


00330353016 


WT 


2008 


10-24 01:56:29 


2008- 


10-24 08:33:59 


3.6 


20.68 



a The exposure time is usually spread over several snapshots (single continuous pointings at the target) during each 
observation. 

b During this observation XRT repeatedly switched between PC and WT modes; since the exposure in WT is short 
(about 39 s), we did not include it in our analysis. 




0.5 1 1.5 2 

Phase 

Figure 4. Time evolution of the 0.5-10 keV light curves folded by means 
of the P—P—P phase-coherent timing solution. Each light curve is covering 
approximately a three days time interval from the BAT trigger, where Tl to 
T7 correspond the the first and last time interval. Light curves are arbitrarily 
shifted along the vertical axis. 



2.2.1 Timing 

The 0.5-10 keV events were used to study the timing properties 
of the pulsar after having corrected the photon arrival times to the 
barycentre of the Solar system with the BARYCORR task (RA=15 fe 
50 m 54! 12, Dec=-54° 18' 24'.' 19 and J2000 was assumed for the 
source position; Israel et al. 2010). We started by inferring an 
accurate period measurement by folding the data f rom the first 
XRT po inting (in WT mode) at the period reported by Cami lo et al.l 
(2007b) and studying the phase evolution within the observation 
by means of a phase-fitting techn ique (details on this technique 
are given in iDalFOsso etal]|2003h . Due to the high variability 
of the pulse shape we decided not to use a pulse template in the 
cross-correlation (which might artificially affect the phase shift), 
but to fit each individual folded light curve with the fundamental 
plus two higher harmonics. In the following we also implicitly 
assumed that the pulsation period (and its derivatives) is a reliable 
estimate of the spin period (and its derivatives), an assumption 
which is usually considered correct for isolated neutron stars. The 
resulting best-fit period (x =1.1 for 2 degree of freedom, hereafter 
d.o. f.) is P = 2.071 3 36(3) s (la al.; epoch 54742.0 MJD; 
see lEsposito et al.ll2008al) . The above period accuracy of 3 /is is 
enough to phase-connect coherently the later Swift pointings car- 
ried out on a daily baseline (see Table[2j. Therefore, the procedure 
was repeated by adding, each time, a further Swift pointing. The 
relative phases were such that the signal phase evolution could be 
followed unambiguously for the whole Swift visibility window (see 
Figure [3). A quadratic term was clearly present in the phase time 
evolution since the very beginning. The resulting phase-coherent 
solution had a best-fit period P — 2.071 336 8(3) s and period 
derivative P = 4.52(4) x 10~ n s s _1 (MJD 54742.0 was used 
as reference epoch; lcr al.). These values are in agreement (within 
~3cr) with those P = 2.071 335(1) s and P = 5.7(5) x 10~ n s 
s~ , reported bv llsraeletaT] ( l2008al) and based on a reduced Swift 
dataset. The fit of the phases with a linear and a quadratic term 
remained statistically acceptable (reduced Xr ~ 0.8 -r 1.2 for 
8-^9 d.o.f.) for the first two weeks since the BAT trigger. However, 
during the latest week the pulse phases increasingly deviate from 
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the extrapolation of the above P—P solution (see Figure[3] central 
panel, empty triangles), resulting in an unacceptable fit (Xr ~ 18 
for 12 d.o.f.). Therefore, we added a higher order component in 
order to take into account for the possible presence of a temporary 
or secular P term. The resulting new phase-coherent solution 
is P = 2.0713401(7) s, P = 3.2(2) x 10 -11 s s -1 , and 
P = 2.0(4) x 10" 17 s s" 2 (MJD 54 742.0 was used as reference 
epoch; la c.L; xl = 2 for 11 d.o.f.), or v = 0.482 779 2(2) Hz, 
v = -7.4(5) x 10~ 12 Hz s"\ and i> = -4.8(9) x 10 -18 Hz 
s -2 . The time residuals with respect to the new timing solution 
are reported in Figure[3] (central panel; empty squares). The 
significance level for the inclusion of the cubic component is 5.1a". 
Moreover, the new timing solution implies a r.m.s. variability of 
only 0.08 s, corresponding to a timing noise level of about 3%, 
usually observed in isolated neutron stars. It is worth emphasizing 
that the second period derivative we found is very unlike to be 
related to a random change of the pulse profiles which is expected 
to introduce only a random distribution of the phase residuals, 
rather then a cubic term. Recent studies on a sample of 366 radio 
pulsars showed that cubic terms on phase residuals are possible 
though smaller by several orders of magnitudes than that detected 
in IE 1547.0-5408 and recorded on time- s cales longer (years) 
than those we are sampling in our dataset iHobbs| [2010). More 
in general the (long-term) timing noise of young radio pulsars 
(t < 10 5 yr) can be exp lained as be ing caused by the recovery from 
previous glitch events (Ho bbs|2010l) . The latter results make unlike 
that the period second derivative is due to random noise. Finally, 
we note that in the above analysis we assumed that the emitting 
region is fixed in time with respect to the o bserver, as suggested by 
studies on other tra nsient magnetars (see IPerna & GotfheHj|2008l ; 
lAlbano et alj l201Cb . Correspondingly, we can also reasonably 
exclude that the cubic term is introduced by a random motion of 
the hot-spot. 
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Figure 5. Swift XRT 2-10 keV observed fluxes (upper panel) as a function 
of the time since the BAT trigger. The two central panels show the evolu- 
tion of the characteristic temperature (kT) and equivalent radius (R) of the 
blackbody model and inferred from the spectral fitting. The lowest panel 
shows the photon index (T) evolution by using a power-law component. 



Based on the above phase coherent solution we also studied 
the pulse shape and pulse fraction evolution. These are changing as 
shown in Figuref3](lower panel) and Figure|4]as a function of time. 
In particular the pulsed fraction (semi-amplitude of the sinusoid 
divided by the mean count rate) has increased smoothly from 20% 
just after the outburst onset, up to almost 50% at the end of the Swift 
visibility window, 22 days after the BAT trigger. The pulse shape 
is variable too, from a multi-peak profile at the beginning of the 
outburst, to a less structured sinusoid during the latest XRT obser- 
vations, though the lower statistics of the latter datasets might have 
hidden a more complex behavior. To assess the significance of the 
observed pulse shape variations as a function of time, we compared 
the folded pairs of contiguous lightcurves reported in Figur e[4]by 
using a two-dimensional K olmogorov-Smimov test dPeacoc kl 19831 ; 
Fasan o~& Franceschinil 1987b . The results show that the probability 
that they do not come from the same underlying distribution is >9a 
(T1-T2), ~2.2a (T2-T3), ~3.3a (T3-T4), ~2.8a (T4-T5), ~3.8a 
(T5-T6), and ~2.9a (T6-T7). 

2.2.2 Spectroscopy 

For the spectral fitting (with XSPEC version 12.4), the data were 
grouped so as to have at least 20 counts per energy bin. The ancil- 
lary response files were generated with XRTMKARF, and they ac- 
count for different extraction regions, vignetting and point-spread 
function corrections. We used the latest available spectral redistri- 
bution matrix (vOl 1) in CALDB. 

Spectral modelling was performed by fitting together all the 



datasets in the 1-10 keV energy range (the high absorption value 
made the data below 1.5 keV almost useless). We first fit the spec- 
tra to simple models: a power law and a blackbody, both modified 
for interstellar absorption. All the parameters were left free to vary, 
except for the absorption column density, that was left free but with 
the request to be the same for all observations. Both models pro- 
vide acceptable fits. The power-law fit (x 2 = 1-04 for 1 292 d.o.f.) 
yields a column density A^h = (6 ± 1) x 10 22 cm -2 and photon 
indices that increase from r ~ 2 to 4 (see Figure[5] lower panel). In 
the data taken immediately after the first BAT trigger, IE 1547.0- 
5408 was found to be in a record high state: the observed flux was 
~7 x 10~ n erg cnT 2 s _1 in the 2-10 keV band, a bout 100 times 
greater than that observed in 2008 July-September (Espos ito et al.l 
2008a ), and ~10 times g reater than the previous highest recorded 
flux faalpern et alj|2008l) . The observed flux declined by 70% in 
three weeks and the fading trend (see Figure|5] upper panel) can be 
described by a power-law with decay index a — —0.171 ± 0.005 
(using the BAT trigger as origin of time). 

From the blackbody fit (x 2 = 1.02 for 1 292 d.o.f.) we derived 
a lower column density of (2.50 ± 0.06) x 10 22 cm -2 . Again there 
is a clear indication of spectral softening during the decay, with the 
blackbody temperature changing from kT ~ 1.3 to 0.7 keV (while 
the radius increased from ~1 to 3 km). 

For those observations for which the statistics was good 
enough (those obtained during the first day since the BAT trigger), 
the Xr was approximately in the 1.2-1.3 range (for 549 d.o.f.), sug- 
gesting a more complex spectrum. In order to check this issue, we 
re-analyzed the first four XRT/PC observations (see TableO fit- 
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ting a BB+PL model. The results are shown in Figure[6]and can be 
summarized as follows: (a) the first spectrum (99 s after the BAT 
trigger) can be totally ascribed to a hard power-law (T = 1.1± 0.2, 
90% c.L), (b) during the second and the third observations (between 
0.05 and 0.2 days after the trigger) the BB component becomes 
dominant (kT = 0.75 ± 0.05 keV in both cases), though a PL 
component is still needed (but not statistically significant) above 
4-6 keV, (c) after one day since the BAT trigger (fourth observa- 
tion) the PL component is not any more detectable in the BB+PL 
model, while in the single PL model the photon index becomes >4 
(mimicking a Wien tail of a BB component). Given the relatively 
low number of counts we did not attempt to fit more sophisticated 
and physical based models. 



3 INTEGRAL OBSERVATION 

Following the flux increase detected by Swift/XKT we triggered 
our IN TEGRAL Target of O pportunity observation program. INTE- 
GRAL dwinkler et atll2003l) observed the source from 2008 Octo- 
ber 8 at 21:35 UT until October 10 at 04 :53 UT for a total exp osure 
time of 95.8 ks. We analyzed the IBIS fobertini et aTll2003r> data, 
the coded mask imager on board INTEGRAL. In particular we se- 
lected the data de rived from its low e nergy (15 keV-1 MeV) detec- 
tor plane ISGRI dLebrun et al.ll2003h . which is the most sensitive 
instrument below 100 keV. 

We analyzed the 48 pointings that composed our observation, 
and derived the corresponding mosaicked images in two energy 
bands: 18-60 keV and 18-100 keV. The source is not detected dur- 
ing our observations, and the 3a upper limits at the source position 
are 0.36 and 0.42 counts s" 1 in the two energy bands respectively. 
This corresponds, assuming a power law spectral shape with photon 
index V — 2.2, to a flux upper limit of 1.8 x 10" 



11 erg cm 2 s 



in the 18-60 keV band, and 2.5 x 10" 11 erg cm -2 s" 1 in the 18- 
100 keV band. Assuming a photon index T = 1.5 the upper limits 
are 2 x 10" 11 erg cm" 2 s" 1 in the 18-60 keV band and 3 x 10" 11 
erg cm" 2 s" 1 in the 18-100 keV energy band. 

In order to search for short bursts, we extracted event lists 
from each pointing. The events were then selected in energy 
(18-60 keV), and extracted from the pixels that were illuminated 
by the source direction for at least 60% of their surface. From these 
event lists, we derived light curves with a binning of 0.05 s. No 
burst was detected above the background level. 



4 DISCUSSION AND CONCLUSIONS 

Thanks to the Swift capability of detecting high energy impulsive 
events with the BAT, and of rapidly pointing to the BAT position 
with the narrow field instrument XRT, we had the chance of study- 
ing the short burst and long-term outbursting properties of one 
of the most interesting objects of the magnetar class, IE 1547.0- 
5408, which entered a renewed phase of activity on 2008 October 
3. In p articular, given the re latively large period derivative of the 
source dCamilo et af] |2007b). we set up an observational strategy 
so that daily pointings were performed during a three week base- 
line (during which the source was observable by Swift). We studied 
the pulsed emission and inferred a rather accurate phase-connected 
timing solution which also includes a second derivative component. 
The spectral evolution was also studied as a function of the decreas- 
ing flux level. 



The connection between the emission of short X-ray bursts 
and the increase of the persistent X-ray emission of AXPs and 
SGRs, has been reliably assessed in the last few years, espe- 
cially thanks to the prompt Swift X-ray follow-up observations 
after the detection of short bursts in CXOU J16471 2-455216, 
SGR 1627-41, SGR 1900+14 and SGR 0501+45 16 ferael et al] 
| 2007al; lEsposito et alj l2008bl : llsrael et alj l2008bl : iHolIand l2008t 
kea et alj2 009). Furthermore, magnetar outbursts have been found 
to be characterized by an intensity-hardness correlation, the harder 
the more intense (and viceversa), as also observed in IE 1547.0- 
5408. Both the connection with the occurrence of the bursts, and 
the hardness intensity correlation point toward a scenario in which 
these events are related to the evolution of a twisted magnetic field, 
which , once implanted by sudden crustal motio ns, necessarily de- 
cays ^Thompson et all2002l ; Heloborodovll2"008h . 

Particularly interesting are the burst energetics recorded by 
BAT. In fact, the 15-150 keV fluence of the bursts reported in Ta- 
bleQJi s in the ~ 1 x 10 39 ergs range (inferred for a distance of 4kpc; 
Tieng oetalJl201Ch . These estimates are conservative lower limits 
to the total energy emitted in the bursts by the source. On the other 
hand, the upper limit on the persistent emission of the source, as 
seen by XRT, is of the order of ~ 4x 10 38 ergs (by integrating the 
exponential flux decay law from the beginning of October 2008 to 
mid- January 2009; see also first panel of Figure[5](, corresponding 
to a ratio between energy released in the bursts and in the persistent 
tail >2. The latter value is in between those observed for SGRs, i.e. 
>10, and those of AXPs where the typical ratio is of the order of 
unity or less. 

The X-ray timing properties of IE 1547.0-5408 in October 
2008 are such that a strong second derivative component, beside 
the usual P and P ones, is needed in order to account for the evo- 
lution of the pulse phases during the about three weeks covered by 
Swift observations. The latter component is P = 2.0(4) x 10" 17 
s s~ 2 or v = —4.8(9) x 10~ 18 Hz s~ 2 , therefore with the same 
sign as that of the first derivative. Second period derivatives have 
been detected only for a few AXPs/SGRs and, in all cases, oc- 
curred injrop_onse_tc^Htctes, s^ of 1RXS J1708- 
4009 i lDall'Osso et al]|2003 : llsrael et alj|2007bl : iDib et alj|2008h 



and IE 1841-045 foib et alJ]2008l) , and/or out bursts, XTE J1810- 
197 dCamilo et alj2007ah and SGR 0501+4516 ( llsrael et al.l2008d : 
iRea et alj200^ . In nearly all cases the second derivatives were rel- 
atively small, while the sign of these components was such that it 
acts in decelerating the spin down rates; typical v values are ap- 
proximately in the 0.05-7-5 xlO" 21 Hzs" 2 range. Only for one 
AXP, namely 1RXS J1708-4009, negative second period deriva- 
tives in the -0.01 -; — 1.3 xlO" 20 Hzs" 2 range has been detected 
just after a glitch foib et al]2008l ; llsrael et al.1l2007bh . 

The timing properties of IE 1547.0-5408 were already known 
to be unusual durin g the first 6 months o f radio monitoring (June 
- November 2007; ICamilo et al.|[2007bh . Even in that case, ob- 
servations were taken after a period of activity and also at that 
time both the first and second frequency derivatives were negative: 
v = -6.5 x 10" 13 Hz s" 1 and v = -1.2 x 10" 19 Hz s" 2 . 
Moreover, during the same time span in which the v (and torque) 
increased, the X-ray flux decreased, a property whi ch is difficult 
to accou nt for with simple assumptions. According to lCamilo et al.1 
(2007b), a possibility is that an extra torque (either than the mag- 
netic spin down) originated from a particle wind that twists the 
magnetic field and, at the same time, transfers some angular mo- 
mentum. In any case, the observed phenomenology calls for a 
mechanism in which X-ray luminosity and spin-down rate are not 
caused by the same mechanism or, at least, not by the same mech- 
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Figure 6. Left panel: Swift XRT spectra and residuals of the first four observations (listed in Table|2j carried out within ~1 day since the X-ray outburst of 
IE 1547.0-5408 (pink filled circles, black, blue and red crosses, from the first to the fourth, respectively). The spectra are fitted by using a blackbody (solid 
lines) plus a power-law (dot-dash lines) model: the hard-to-soft evolution is evident (see text for details). Right panel: same as before but for the density flux 
in order to emphasize the spectral evolution above ~4-5 keV. 



anism and also at the same time. Even more dramatic is the case of 
October 2008 reported in this paper, with IE 1547.0-5408 show- 
ing a v which is a factor of about 50 larger than in 2007, while 
the X-ray flux is a factor of 4-30 higher than during 2007. It is 
very likely that the v components related with glitch-only events 
have a different origin than those related to outburst events. For 
another transient AXP, XTE J1810-197, a v component has been 
det ected during an outb urst, although after two years since the on- 
set dCamilo et alj|20075) . Correspondingly, a straightforward com- 
parison with the IE 1547.0-5408 properties is not feasible (similar 
arguments also hold for the radio properties of the two sources). In 
the case of SGR 050 1+45 16, the second derivative component was 
inferred during the first two months since the outburst onset and 
has a sign opposite to that of IE 1547.0-5408. 

A plausible scenario is that the dissimilar behavior seen in 
the various sources is related to the difference between the out- 
burst P ou tb and the pulsar secular one P SO c- Regardless of the 
mechanism responsible for the outburst and for the net timing ef- 
fects at the outburst onset, it is reasonable to imagine the pulsar as 
a system which tends to recover its (almost) initial ti ming prop- 
erties, i. e. same pre-outbu rst and post-outburst P scc dAlpar et aU 
1984allb1: lLmketal.lll992h . Correspondingly, a negative or posi- 
tive second period derivatives might be ascribed to whether the 
ioutb over Psec ratio is larger or smaller than unity. Unfortunately, 
for all the four magnetars for which an outburst was detected in 
the latest few years (XTE J1810-197, CXOU J164710.2-455216, 
SGR 0501+4516 and IE 1547.0-5408) the P scc is unknown. In the 
case of IE 1547.0-5408 the positive second derivative detected in 
2007 and 2008 would imply a P out b smaller than P scc '- the lat- 
ter being the value to which the source is approaching asymptoti- 
cally. Within the classical scenario of isolated neutron stars, P sec 
is the spin-down value to be used when inferring the pulsar dipo- 
lar magnetic field strength. Correspondingly, the inferred value of 
Bd ~ 2.7 x 10 14 Gauss should be considered as a lower limit on 
the intrinsic dipolar magnetic field strength. 

Taken at face values, the increase of P during the outburst de- 
cay appears difficult to reconcile with the predictions of the twisted 
magnetosphere model. In such model, in fact, the spin-down rate 



is dictated by the radial variation of the external field. A globally 
twisted dipolar field has radial dependence ~ R~ 2 ~ p , where the 
index p decreases with increasing twist angle (0 ^ p ^ 1). If, as 
it seems natural to assume, the twist decays as the source returns 
to quiescence, one expects a decrease in the spin- down rate, i.e. a 
negative P. As pointed out by Belob orodov (20083, however, in re- 
alistic situations the twist is likely to be confined to a rather narrow 
bundle of current-carrying field lines. In this case, and if the twist 
is suddenly implanted but its strength is moderate (twist angle less 
than ~1 rad), the twist may still grow for a while in spite of the 
luminosity released by dissipation monotonically decreasing. Only 
after the twist angle has reached its maximum value, it will start 
to decrease, together with the torque. This implies that the spin- 
down rate may be non-monotonic during the outburst decay, with a 
first phase characterized by a positive P followed by a "recovery" 
during which P < and P goes back to its secular value. The typ- 
ical duration of the increasing spin-down stag e is difficult to quan - 
tify, but it can last ~ months (see Eq. 55 in lBeloborodovll2008l). 
Indeed , this effect has been already proposed by iBeloborodovl 
(2008) to explain the anti-correlation between X-ray lumino sity 
and torque observed in IE 1048.1-5937 dGavriil & KaspHl2004l) . In 
IE 1547.0-5408 P has been observed to increase monotonically 
during summer/fall 2007 ( during the decay of an outburst the peak 
of which has been missedi lCamilo et alj2008l). The value of ii cor- 
responding to the latest radio observation of ICamilo et al.l j2008t) 
is ~ —7.35 x 10~ 12 s~ 2 , and that in October 3, at the beginning 
of the X-ray outburst, ~ -(7.4 ± 0.4) x 10~ 12 s~ 2 . A continu- 
ous decrease in v seems therefore hardly compatible with observa- 
tions, and, since a constant v has never been observed in AXPs, it 
may well be that v started increasing at some epoch between 2007 
November 22 (latest radio observations) and 2008 October 3 (X-ray 
outburst onset). This behavior would not be in contradiction with 
the untwisting magnetosphere model mentioned above. Finally, it is 
worth mentioning the possibility that the unusual P behaviour after 
the October 2008 outburst might be one of the reasons for the ex- 
tremely intense outburst displayed by IE 1547.0-5408 in January 
2009. In other words, the positive P would be acted in accelerating 
the onset of stresses on the NS surface. 
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The pulsed fraction of IE 1547.0-5408 clearly increases dur- 
ing October 2008, a behavior already observed during the initial 
phases of AXP outbursts, such as the case of CXOU J 1647 10.2- 
455216 where the quiescent pulsed fraction (~80%) dropped to 
~10% at the outburst ons et and started incr easing thereafter to- 
wards the quiescent value dlsrael et alj|2007^ similar results have 
been recently obtained for SGR 0501+45 16: lRea et al.120091) . How- 
ever, in the case of IE 1547.0-5408 the pulsed fraction at the out- 
burst onset (~20%) is already larger than the quiescent value (as 
seen by XMM-Newton; ~7%) and it reached ~50% at the end of 
the October 2008 Swift monitoring window. This finding, together 
with the very variable pulse profile as a function of time, suggests 
that the emission and/or magnetospheric geometry is complex and 
variable during the initial phases (on month timescales) of the out- 
burst, and likely not easily related to the geometry in quiescence. 



Pulsations in the X-ray flux may be a consequence of emission 
from a limited area of the star surface, of non-isotropic magneto- 
spheric emission, or a combination of both effects. Again, if the 
star surface is heated by the returning currents which appear when 
a magnetospheric twist is implanted, even considering the simple 
evolution of the twist dictated by ohmic dissipation, an increase of 
the pulsed fraction in time is expected because, as time elapses, the 
bundle of current-carrying field lines shrinks and the heated area 
decreases. On the other hand, an increase in both pulsed fraction 
and emitting area, which is not excluded by the present data, would 
point toward a complete disassociation between the hot spot and 
the timing properties, with the latter perhaps of magnetospheric ori- 
gin. The XRT statistics is not enough to unambiguously assess the 
evolution of the emitting area, nor even the presence of two dis- 
tinct spectral components (BB and PL) detected by XMM-Newton 
in 2006 and 2007, although there are several hints suggesting that 
also in the spectra observed during October 2008 two components 
are needed. 



Overall, the spectral properties of IE 1547.0-5408 are also 
very similar to those of other transient magnetars, with an evident 
decay of the flux and a softening of the emitted photons (at least in 
the 1-10 keV range) as the outburst evolves in time. The spectral 
softening is apparent while using both the BB and the PL spectral 
models. Referring, for concreteness, to the power-law model, the 
photon index switches from ~ 2 at the outburst onset to ~ 4 af- 
ter ~ 10 d. A softening is expected if the spectrum originates in a 
globally twisted magnetosphere by resonant up-scattering. As the 
twist angle decreases the charge density decreases making scatter- 
ing less efficient. This implies that the high-energy tail filled up 
by up-scattered thermal photons becomes less and less populated, 
producing a softening of the spectrum. This has been verified by 
detailed Monte C arlo simulations in the case of a globally twisted 
dipolar field (e.g. iNobili et al.ll2008l) . Although no detailed model 
has been presented as yet, it seems likely that the same picture holds 
if a localized twist is considered instead. As it was mentioned ear- 
lier, the twist angle may still increase for a while following the 
onset of the outburst, which would b e consistent with th e negative 
P, while its spatial extent decreases dBeloborodovll2008r) . This lat- 
ter effect may compensate for the first one and be the main effect 
responsible for the spectral evolution: in particular, if the region 
in which photons can be efficiently scattered by the currents is re- 
duced, a flux softening and decrease in luminosity are expected. 
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